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Abstract 



We have searched a sample of 9.6 million BB events for the flavor-changing 
neutral current decays B ^ K£+£- and B K* {892)1+ £~. We subject the 
latter decay to the requirement that the dilepton mass mu exceed 0.5 GeV. 
There is no indication of a signal. We obtain the 90% confidence level upper 
limits I3{B Ki+e-) < 1.7 x 10"^ and 13{B K*(892)£+£-)^,,>o.5GeV < 
3.3 X 10^^. We also obtain an upper limit on the weighted average 0.6513{B 
Ki+e~-) + 0.35B{B K*{892)i+£-)mu>o.5GeV < 1-5 x 10"^ The weighted- 
average limit is only 50% above the Standard Model prediction. 



1 



Typeset using REVT^ 



S. Anderson/ V. V. Frolov/ Y. Kubota/ S. J. Lee/ R. Poling/ A. Smith/ 
C. J. Stepaniak/ J. Urheim/ S. Ahmed/ M. S. Alam/ S. B. Athar/ L. Jian/ L. Ling/ 
M. Saleem/ S. Timm/ F. Wappler/ A. Anastassov/ E. Eckhart/ K. K. Gan/ C. Gwon/ 
T. Hart/ K. Honscheid/ D. Hufnagel/ H. Kagan/ R. Kass/ T. K. Pedlar/ J. B. Thayer/ 
E. von Toerne/ M. M. Zoeller/ S. J. Richichi/ H. Severini/ P. Skubic/ A. Undrus/ 
V. Savinov/ S. Chen/ J. W. Hinson/ J. Lee/ D. H. Miller/ E. I. Shibata/ 
I. P. J. Shipsey/ V. Pavlunin/ D. Cronin-Hennessy/ Y. Kwon/'[| A.L. Lyon/ W. Park/ 
E. H. Thorndike/ T. E. Coan/ Y. S. Gao/ Y. Maravin/ 1. Narsky/ R. Stroynowski/ 

J. Ye/ T. Wlodek/ M. Artuso/ K. Benslama/ C. Boulahouache/ K. Bukin/ 
E. Dambasuren/ G. Majumder/ R. Mountain/ T. Skwarnicki/ S. Stone/ J.C. Wang/ 

A. Wolf/ S. Kopp/° M. Kostin/° A. H. Mahmood/^ S. E. Csorna/^ I. Danko/^ 
K. W. McLean/2 Z. Xu/^ R. Godang/^ G. Bonvicini/^ D. Cinabro/^ M. Dubrovin/^ 
S. McGee/^ A. Bornheim/^ G. Eigen/^'f] E. Lipeles/^ S. P. Pappas/^ A. Shapiro/^ 
W. M. Sun/5 A. J. Weinstein/5 D. E. Jaffe/^ R. Mahapatra/*^ G. Masek/^ H. P. Paar/*^ 

A. Eppich/^ R. J. Morrison/'^ R. A. Briere/^ G. P. Chen/*^ T. Ferguson/^ H. Vogel/^ 
J. P. Alexander/^ C. Bebek/^ B. E. Berger/^ K. Berkelman/^ F. Blanc/^ V. Boisvert/^ 
D. G. Cassel/9 P. S. Drell/^ J. E. Duboscq/^ K. M. Ecklund/^ R. Ehrlich/^ 
P. Gaidarev/9 L. Gibbons/^ B. Gittelman/^ S. W. Gray/^ D. L. Hartill/^ 
B. K. Heltsley/^ L. Hsu/^ C. D. Jones/^ J. Kandaswamy/^ D. L. Kreinick/^ M. Lohner/^ 

A. Magerkurth/9 H. Mahlke-Kriiger/^ T. O. Meyer/^ N. B. Mistry/^ E. Nordberg/^ 
M. Palmer/^ J. R. Patterson/^ D. Peterson/^ D. Riley/^ A. Romano/^ H. Schwarthoff/^ 
J. G. Thayer/^ D. Urner/^ B. Valant-Spaight/^ G. Viehhauser/^ A. Warburton/^ 
P. Avery/° C. Prescott/° A. I. Rubiera/° H. Stoeck/° J. Yelton/° G. Brandenburg/^ 
A. Ershov/i D. Y.-J. Kim/^ R. Wilson/^ B. I. Eisenstein/^ J. Ernst/^ G. E. Gladding/^ 
G. D. Gollin/2 R. M. Hans/^ E. Johnson/^ I. Karliner/^ M. A. Marsh/^ C. Plager/^ 
C. Sedlack/2 M. Selen/^ J. J. Thaler/^ J. WiUiams/^ K. W. Edwards/^ A. J. Sadoff/^ 
R. Ammar/5 A. Bean/^ D. Besson/^ and X. Zhao^^ 

^University of Minnesota, Minneapolis, Minnesota 55455 
^State University of New York at Albany, Albany, New York 12222 
^Ohio State University, Columbus, Ohio 43210 
^University of Oklahoma, Norman, Oklahoma 73019 
^University of Pittsburgh, Pittsburgh, Pennsylvania 15260 
^Purdue University, West Lafayette, Indiana 47907 
''University of Rochester, Rochester, New York 14627 
^Southern Methodist University, Dallas, Texas 75275 
^Syracuse University, Syracuse, New York 13244 
^'^University of Texas, Austin, Texas 78712 
^^University of Texas - Pan American, Edinburg, Texas 78539 
^^Vanderbilt University, Nashville, Tennessee 37235 



*Permanent address: Yonsei University, Seoul 120-749, Korea. 
^Permanent address: University of Bergen, 5007 Bergen, Norway. 



2 



i^Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061 
^"^Wayne State University, Detroit, Michigan 48202 
^^Cahfornia Institute of Technology, Pasadena, California 91125 
^^University of California, San Diego, La Jolla, California 92093 

^^University of California, Santa Barbara, California 93106 
^^Carnegie Mellon University, Pittsburgh, Pennsylvania 15213 
Cornell University, Ithaca, New York 14853 
^'^Univcrsity of Florida, Gainesville, Florida 32611 
^^Harvard University, Cambridge, Massachusetts 02138 
^^University of Illinois, Urbana-Champaign, Illinois 61801 
^^Carleton University, Ottawa, Ontario, Canada KIS 5B6 
and the Institute of Particle Physics, Canada 
^^Ithaca College, Ithaca, New York 14850 
^^University of Kansas, Lawrence, Kansas 66045 



3 



The flavor-changing neutral current (FCNC) decay b si~^i~ is sensitive to physics beyond the 
Standard Model Q, and, like the radiative penguin decay b — s- ,57, is more amenable to calculation 
than purely hadronic FCNC decays. The decay b — > si~^i~ depends on the magnitude and sign of 
the three Wilson coefficients C7, Cg, and Cio in the effective Hamiltonian, while 6 — > 57 depends 
only on the magnitude of C7. These three Wilson coefficients are likely places for New Physics 
to appear, as they come from loop and box diagrams. Upper limits on the branching fraction for 
b — s- si~^£~ thus place constraints on New Physics, while observation of 6 — s- s£'^£~ at a rate in 
excess of that predicted by the Standard Model would provide evidence for New Physics. Just as we 
first observed 6 — > 57 through its exclusive decay B — > (892)7 Q, and only later in an inclusive 
fashion [g, so we search for b si+i~ through the decays B Kl+l' and B K*{892)£+e-. 
Existing published upper limits come from CDF [Q], and from earlier work by us Here we 
present improved upper limits. 

The B K*£~^£~ decay rateQ peaks at low dilepton mass mu, due to the photon pole from 
B K*^mrtuaii Ivirtuai £'^ £~ ■ Because the decay B — K*^ is already well studied, and IC7I 
thus reasonably well known, we require B — > candidates to have dilepton mass above 0.5 

GeV, to reduce the contribution from the virtual photon diagram and thus the dependence on \Cj\. 

The data used in this analysis were taken with the CLEO detector [g at the Cornell Electron 
Storage Ring (CESR), a symmetric e'^e~ collider operating in the T(4S) resonance region. The 
data sample consists of 9.2 fb~^ at the resonance, corresponding to 9.6 million BB events, and 
4.5 fb~^ at a center-of-mass energy 60 MeV below the resonance. The sample below the resonance 
provides information on the background from continuum processes e'^e~ qq, q = u,d, s, c. 

We search for B — > K^*^£^£~ both in the ^"''z^" and e~^e~ modes, for B — > K£^£~ in both 
the and modes, and for B —>■ K*£^£~ in the K*^ — > K~^Tr~ and K^tt^ modes and in 
the K*^ — > K^TT^ and K^tt^ modes, a total of 12 distinct final states. is detected via the 
tt'^tt decay chain. 

There are three main sources of background: B K(*^('\ V^'^ ^ and other B 

decays; continuum processes with two apparent leptons (either real leptons or hadrons misidentified 
as leptons); and BB decays other than B ip^'^X, with two apparent leptons. We suppress B — > 
^{*)^{') events with vetoes on the dilepton mass around ^ and ip' , in particular 2.80 < niee < 3.23 
GeV, 2.90 < nif^f, < 3.20 GeV, 3.51 < niee < 3.77 GeV, and 3.55 < m^^ < 3.74 GeV. These very 
wide cuts are needed because of the low-side radiative tail from internal and external bremsstrahlung 
from ip^'^ — > e~^e~ , and to a lesser extent the low-side radiative tail from internal bremsstrahlung 
from the li'^ ii~ decay. 

We suppress the background from BB semileptonic decays with a cut on the event missing 
energy, Emiss, since events with leptons from semileptonic B or D decay contain neutrinos. Dis- 
tributions in Emiss for Monte Carlo samples of signal and background are shown in Fig. ||. We 
suppress continuum events with a cut on a Fisher discriminant, a linear combination of R2 (the 
ratio of second and zeroth Fox- Wolfram moments Q of the event), cos 6tt {Ott the angle between the 
thrust axis of the candidate B and the thrust axis of the rest of the event), S (the sphericity), and 
cos 9b {Ob the production angle of the candidate B, relative to the beam direction). In particular, 
= ^2 + 0.1171 cos Qui + 0.779(1 - S) + 0.104| cos Ob]- The coefficients of all terms but R2 were 



*Throughout this Letter, the symbol K* means i^*(892). 
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determined by the standard Fisher discriminant procedure The relative weight given to R2 
was determined visually, from a scatter plot of R2 vs. the Fisher discriminant from the other three 
variables. Distributions in for Monte Carlo samples of signal and background are shown in Fig. |l|. 

For those decay modes involving a charged kaon, we use specific ionization {dE/dX) and time- 
of-flight information to identify the kaon, cutting loosely (3 standard deviations) if those variables 
deviate from the mean for kaons in the direction away from the mean for pions, and harder (by a 
variable number, denoted klD'^'^^, of standard deviations) if they deviate on the side towards the 
pions. 

All cuts have been determined from Monte Carlo samples: continuum events, BB events with 
no signal, and BB events with a B ^ K^*^£^(.~ signal. We optimized cuts on !F, E^igg, and, where 
appropriate, kaon identification simultaneously. We found that the optimum curve in efficiency vs. 
background space was traced out if we required that cuts on J^, Emissi ^-nd kID were tightened 
or loosened together. In particular, we found that the optimum curve was well described by 
E^l^ = 1.0 + 6.67 X (J^^-^"* - 0.8) GeV, A:IL>^"* = 3.0 - 0.66 x (3.0 - E^lJ. With this formulation, 
we have a single cut variable, .F'^"*, to optimize. 

We optimized the cuts to obtain either the best upper limit, assuming no signal, or to see the 
smallest possible signal. These two different optimization procedures led to similar cuts, and we 
took the average. In the optimization we allowed the value of J^^'^^ to vary from decay mode to 
mode. The final cuts on T are shown in Table |. The corresponding cuts on Emiss and kID can 
be obtained from the expressions given above. 

Our final discrimination between signal and background comes from the B reconstruction 
variables conventionally used for decays from the T(4S'), beam-constrained mass M^and = 

\J ^beam ~ ^cand IS^E = Ecand - Eheam- Our resolution in Mcand IS 2.5 MeV, and in A£', 
20 MeV. We define a signal box in Mcand — space, ± 6.5 MeV around the B mass by ± 60 
MeV (/i^/x~) or it 70 MeV {e^e~). The signal box for e^e~ events is shifted off zero by 5 MeV in 
Ai?, because radiation losses cause B i^^*^e^e~ signal events to peak at -5 MeV in /S.E rather 
than at zero. 

Background from B ip^'^X is estimated from Monte Carlo simulation. Background from 
other B decay processes and from continuum processes is determined using a large sideband region 
in Mcand - A£' space: 5.20 < Mcand < 5.29 GeV, \AE\ < 0.25 GeV, but excluding the signal 
box. From Monte Carlo simulation, we found that the ratio of events in the signal region to events 
in the sideband region is 0.024 for BB background events other than B — > ip^'')X, and 0.027 for 
continuum background, in both cases smaller than the ratio of areas, 0.038, because backgrounds 
fall off as the candidate mass approaches the beam energy. Recognizing that this ratio must be 
larger for continuum background than for BB background because the more jet-like continuum 
events will not fall off as rapidly as candidate mass approaches beam energy, and recognizing that 
using a lower background in an upper limit computation gives a more conservative answer, we take 
the continuum scaling factor to be equal to the BB scaling factor, rather than using 0.027. 

The e~^e~ events suffer a degradation in resolution due to internal and external bremsstrahlung 
from the electrons. We partially recover that resolution by adding to each electron energy the 
energy of those photons found nearby in angle. This procedure improves the '0 veto, and resolution 
in Emiss, Mcand, and AE. 

The number of events that satisfy all cuts and land in the signal box is given, for each mode, in 
Table ||, along with the background estimate. We find 3 — > K£~^£~ candidates, with an expected 
background of 2.0; we find 4 — > K*£~^i~ candidates, with an expected background of 3.8. Thus 
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there are a total of 7 events with an expected background of 5.8. The probabihty that a true mean 
of 5.8 win fluctuate up to 7 or more events is 36%. Thus, there is no indication of signal. A scatter 
plot of Mcand VS. AE for events passing all other cuts and landing in the signal or sideband region 
is shown in Fig. ||. Again, no indication of a signal. We obtain upper limits. 

We calculate upper limits, at 90% confidence level, taking backgrounds into account To 
allow for the uncertainty in the background estimate, we use a value for the background which is 
reduced below our actual estimate. For individual modes, we use half the estimated background. 
For B — i- Ki~^i~ and B — > K*£^i~ totals, we use the estimated background minus 1.28 standard 
deviations of its combined statistical and systematic error. (The factor 1.28 gives a 90% confidence 
level lower estimate of the background, assuming its uncertainty has a Gaussian distribution.) 
Results are given in Table |. 

We gain experimental sensitivity to the underlying b — > si~^i~ interaction by calculating a 
weighted average over the two decay modes studied. To account for the difference in the exper- 
imental precision for each mode, we weight them by our relative efficiencies, that is, we com- 
pute an upper limit on the sum over all the individual sub-modes. This gives an upper limit on 
OmB{B K£+£-) + 0.35B{B K* £+£-), where the coefficients 0.65 and 0.35 are the relative 
efficiencies we have for the two modes. 

We use Monte Carlo simulation to determine the efficiency for detecting the signal modes. 
The decays B K£~^£~ and B — > K*£~^£~ are generated using the model of Ali et al. [|l|. The 
helicity of the K* is taken into account. Final-state radiation is included, using the CERNlib 
subroutine Photos [jlO|. We have also generated decays with the two extreme variations that Ali 
et al. suggest for their model, and with several other models |11,|]. We find the model-to-model 
variation in efficiency to be small, with a relative r.m.s. variation of ±3%. 

To check our procedures, we have looked for the decays B J/ifjK^*'^ rather than vetoing 
them. We compare the branching fractions obtained with those from prior CLEO measurements. 
There is good agreement - differences are at or below the one-standard-deviation level. 

Systematic errors are of two varieties - those on the estimate of signal detection efficiencies, 
and those on the estimate of backgrounds. The contributors to the former are lepton identification 
uncertainties (contributing ±5%, relative, in the efficiency), missing-energy-simulation uncertainties 
(±3.5%), and simulation uncertainties for Mcand-, AE, dE/dX, time of flight, and (±3.0%), giving 
a ±7% relative uncertainty in the overall efficiency. To this we add in quadrature ±3% for the model 
dependence of the efficiency, discussed earlier. The contributors to the background uncertainties 
are the modelling of 5 — > ip^'^X (±10%), and uncertainties in the scale factors from sideband region 
to signal region in Mcand — AE space. We assign a systematic error to the BB background scale 
factor by determining it with different methods, obtaining 0.024 ± 0.004. Recognizing that the 
scale factor for continuum should be larger than that for BB, we conservatively set it equal to the 
BB scale factor, with the same (correlated) systematic error. The errors shown on the backgrounds 
in Table | include statistical errors and the systematic errors just described. 

There is no universally agreed-upon procedure for including systematic errors in upper-limit 
estimates. We conservatively reduce the background by 1.28 standard deviations, and decrease the 
efficiency by 1.28 standard deviations. In this way we obtain our final results: 
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B{B Ktr) < 1.7 X 10"^ , 
B{B ^ i^*(892)^+r)™„>o.5GeV < 3.3 x 10"^ , and 

0.65S(5 ^ Kl+r) + 0.35i3(S ^ K*(892)£+r )m,,>o.5GeV < 1-5 x 10"*^ , 

all at 90% confidence level. These results are significant improvements over previously published 
limits 

The Standard Model values for these branching fractions, as given by Ali et al. |jl], are 0.6 x 10~^ 
for B{B Ki+i-) and 1.8 x 10"^ for B{B K*l+r)mu>o.bGeV, and thus 1.0 x 10^^ for the 0.65 
/ 0.35 weighted average. The limit on the branching fraction for B — > Ki~^i~ is therefore about 
three times its Standard Model prediction, the limit on the branching fraction for B K*i~^i~, 
subject to the requirement that ma > 0.5 GeV, is about twice its Standard Model prediction, and 
the 0.65 / 0.35 weighted average is only 50% larger than its Standard Model prediction. 

In summary, we have searched for the decays B K£^£^ and B — > i^*(892)^~'"^~. We find 
no indication of a signal, and obtain upper limits on the branching fractions. These limits are 
consistent with Standard Model predictions, but not far above them. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent luminosity 
and running conditions. This work was supported by the National Science Foundation, the U.S. 
Department of Energy, the Research Corporation, the Natural Sciences and Engineering Research 
Council of Canada, the Texas Advanced Research Program, and the Basic Science program of the 
Korea Research Foundation. 
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FIG. 1. Distributions in Emiss (upper) and (lower) for Monte Carlo samples of signal events 
(solid) , BB background events (dotted) , and continuum background events (dashed) . The vertical 
scale is arbitrary. 
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FIG. 2. Scatter plot of Mcand vs. AE for ee events (open circles) and fifi events (solid squares) 
passing all other cuts, for data on the T{4S) resonance. The smaller box (solid) is the ^~ signal 
box, while the larger box (dashes), shifted 5 MeV toward negative Ai?, is the e+e~ signal box. 
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mode 


j:cut 


observed 
events 


background 


efficiency 


B X 10^ 
upper Lim. 




0.938 


1 


0.10 


0.053 


7.6 




0.925 





0.21 


0.041 


7.8 


K+e+e- 


0.925 


1 


0.95 


0.165 


2.3 


K+irir 


0.850 


1 


0.7-i 


0.111 


3.1 


Ki+i- 




3 


1.99 ±0.35 


0.370 


1.49 




0.925 





0.35 


0.019 


12.8 




0.900 





0.27 


0.015 


15.6 




0.800 


3 


0.27 


0.015 


46.0 




0.750 





0.49 


0.008 


29.3 


— + — 


0.925 


1 


0.97 


0.071 


5.0 




0.875 





1.24 


0.052 


4.6 




0.900 





0.11 


0.007 


35.8 




0.750 





0.10 


0.002 


117.3 


K*e+e- 




4 


3.80±0.57 


0.188 


2.94 


Sum 
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5.79±0.83 


0.558 


1.35 



TABLE I. Value of T'^'^^ , number of events observed in signal window, the expected background, 
efficiency, and upper limit on the branching fraction, for the 12 modes, and for B Kt^t~ ^ 
B K*e+i-, and 0.65^(5 Ki+i-) + 0.35B{B K*i+i~). Note, the upper limits given in this 
table have been corrected for systematic error in background, but not corrected for the systematic 
error in efficiency. 
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